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I derive a loop representation for the canonical and grand-canonical partition functions for an 
interacting four-component Fermi gas in one spatial dimension and an arbitrary external potential. 
The representation is free of the "sign problem" irrespective of population imbalance, mass imbal- 
ance, and to a degree, sign of the interaction strength. This property is in sharp contrast with the 
analogous three-dimensional two-component interacting Fermi gas, which exhibits a sign problem 
in the case of unequal masses, chemical potentials, and repulsive interactions. The one-dimensional 
system is believed to exhibit many phenomena in common with its three-dimensional counterpart, 
including an analog of the BCS-BEC crossover, and nonperturbative universal few- and many-body 
physics at scattering lengths much larger than the range of interaction, making the theory an inter- 
esting candidate for numerical study. Positivity of the probability measure for the partition function 
allows for a mean-field treatment of the model; here, I present such an analysis for the interacting 
Fermi gas in the SU(4) (unpolarized, mass-symmetric) limit, and demonstrate that there exists a 
phase in which a continuum limit may be defined. 
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I. INTRODUCTION 



In recent years, nonrelativistic Fermi gases have become the focus of intense study, fueled primarily by significant 
advances in ultra-cold atom experiments [l| . Of primary interest has been the BCS-BEC crossover regime for a dilute 
two-component Fermi gas with an attractive short-range two-body interaction in three spatial dimensions. In the 
limit that the scattering length is much larger than the range of interactions, this system, often referred to as a 
"unitary Fermi gas," becomes universal in the sense that the long distance physics is independent of the details of the 
inter-particle potential. In the unitary limit, the Fermi gas is characterized by a single dimensionful scale, namely the 
density of the system, and as a consequence, the energy of the interacting gas is proportional to the free-gas energy, 
and the pairing gap for fermion excitations is proportional to the Fermi energy. This conformal and scale invariant 
system is not only experimentally realizable using ultra-cold atoms tuned to a Feshbach resonance, but also provides 
an idealized theoretical description of dilute nuclear matter . 

The unitary Fermi gas is a challenging system to explore theoretically due to the strongly coupled nature of the 
interaction, as well as a lack of scales in the problem. Although qualitatively much can be said about unitary fcrmions, 
a nonperturbative treatment is essential for achieving a complete quantitative understanding of the system. Through 
the use of numerical techniques, considerable headway has been made in exploring the properties of the few- and 
many-body system both at zero and at finite temperature. Monte Carlo techniques have been particularly successful 
in studies of the unpolarized unitary Fermi gas, in which each spin component has the same mass (for recent studies, 
see e.g., [H-Q), and the slightly polarized gas, in which a single fermion has been removed [Tf. 

Introducing a population imbalance in the system by means of unequal chemical potentials for each species in 
the grand-canonical ensemble, a mass imbalance, or repulsive interaction results in a sign or signal/noise problem, 
making numerical calculations which relv on Monte Carlo simulations a formidable challenge. Despite the considerable 
attention such systems have received Q, from a theoretical standpoint less is known quantitatively about the phase 
diagram in each of these scenarios. Recently, progress has been made from the experimental side, however The 
mass-imbalanced unitary Fermi gas is particularly interesting because the scale-invariance of the system is expected 
to break down to a discrete scale invariance at a critical mass ratio, leading to Efhmov physics (ill. [T2j|. 

The realization of optical lattices in recent years make experimental exploration of Fermi gases in lower and mixed 
spatial dimensions a feasible and attractive alternative . Theoretical investigations suggest that such systems can 
also exhibit scale invariance and universal properties, much in the same way as their two-component, three-dimensional 
Fermi gas counterpart [l^ . One specific example, which will be the focus of this paper, is the four-component Fermi 
gas in one spatial dimension with an attractive short-range four-body interaction. It was shown that this system is 
also stable and universal in the sense that its properties are independent of the confining potential when tuned to a 
four-body resonance [l5| . 

The relevant dimensionful scales which describe the one-dimensional Fermi gas are the Fermi momentum p^? = 7m/4, 
where n is the number density of the system, the Fermi energy Ep = pjp/ (2m), where m is the mass of the fcrmions, 
and the scattering length a, which characterizes the long-distance properties of the four-body interaction. As with the 
universal Fermi gas in three dimensions, the properties of this model are smoothly controlled by a single dimensionless 
parameter given by 77 = {ppa)~^ [Tsj . In the regime 77 <^ —1, the system forms a one-dimensional analog of a BCS 
superfluid, with three gaps for fermion excitations, each of which are exponentially small in 77. On the other hand, 
in the regime 77 ^ 1, the many-body system forms a dilute Bose-Einstein condensate of deeply bound tetramers with 
largely gapped fermion excitations proportional to l/(r7ia^). Each of these regimes may be understood within the 
framework of effective field theories; in the former case, the low-energy constants entering the effective theory can be 
matched onto the underlying microscopic theory perturbatively in |7y|~^, but in the latter case, the matching must be 
done nonperturbatively by solving the eight-body problem [l5|. 




In the cross-over region corresponding to I77I <^ 1, there remains only one dimensionful parameter to describe 
the system, namely the density. Consequently, by dimensional analysis, the ground state energy density must be 
proportional to the free-gas energy density, Spc = nEp/3. Similarly, the excitation gaps are in proportion to the 
Fermi energy. As is the case for the two-component three-dimensional Fermi gas, the four-component one-dimensional 
Fermi gas satisfies universal relations (i.e.. Tan relations) involving a recently introduced parameter, known as the 
"contact density" [l6l - [20j . which is related to the fall-off of the tail of the momentum distribution of the Fermi gas. 
In addition, as a result of an operator-state correspondence [2lj, an exact relationship exists between the energies of 
the Fermi gas at infinite scattering length when confined to an external harmonic trap, and the scaling dimension of 
primary operators in free-space, thus motivating nonperturbative studies of the trapped few-body system. 

Because the one-dimensional Fermi gas is nonperturbative in the cross-over regime, numerical simulations are crucial 
for obtaining a quantitative understanding of the few- and many-body physics of this model. In this paper, I show that 
an effective theory description of the interacting four-component Fermi gas in one spatial dimension can be formulated 
on a lattice in a way that avoids the sign problem. Surprisingly, the sign problem remains absent in the new path- 
integral representation for unequal masses and chemical potentials, in the presence of an arbitrary external potential. 
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and even to an extent, irrespective of the sign of the four-body interaction strength^ The canonical and grand- 
canonical partition functions are formulated in terms of a path-integral over self-avoiding non-intersecting fermion 
loops, and are suitable for numerical study using a simple loop-updating scheme, or more sophisticated algorithms. I 
derive expressions for both of these partition functions, as well as for expectation values of the kinetic, potential and 
interaction energies, charges and their associated susceptibilities, and iV-point correlation functions involving local 
density operators. Finally, I explore the number density phase diagram of the lattice theory at zero temperature, and 
demonstrate that, at least within mean-field theory, a continuum limit exists. 



II. LATTICE THEORY 



In the universal regime, the thermodynamics of a dilute one-dimensional four-component Fermi gas may be reliably 
described by the continuum effective field theory: 



ScontU-i) = / (i^X 



2m 



(1) 



defined in Euclidean spacetime. Here, ^ = {4'a,'4'bi'4'c,'4'd) is a four-component Grassmann-valued spinor, and m = 
diag (ma, rrib, TTic, md) and ^ = diag (/ia, /i;,, /Xc, /z^) are diagonal mass and chemical potential matrices, respectively. 
Fermion interactions are achieved via a four-fermion contact potential, with strength g tuned to reproduce the physical 
four-body scattering length. I discretize Eq.^on an x Ns lattice A, imposing anti-periodic boundary conditions in 
time and open boundary conditions in space 0. The lattice sites are labeled by the integer coordinate pair n — {ur, Us) 
with rir S [0, A^T- — 1] and Us G [0, A^ — 1]. The physical length of the lattice is given by L = bsNg and the inverse 
temperature by /3 = brNr, where bs and br are the spatial and temporal lattice spacings, respectively. Following 
the approach of [l^, the lattice-discretized action is obtained by replacing the continuous fields by their discrete 
counterpart: ip{x) — > 'ipn, the integral in Eq. [1] with a discrete sum over lattice sites: J <Px ^t&s X^neA' ^^"^ 
identifying the terms in the action as follows: 



1 

Vr 
1 



(2) 



where are unit vectors in the time (j = 0) and space (j = 1) directions of the lattice. The continuum limits in 
space and time correspond to taking all scales with dimensions of length much larger than bg and all scales with 
dimensions of time much larger than b^ ■ For simplicity, I will work in "lattice units" where the temporal and spatial 
lattice spacings (along with h) are set to unity, and where all other dimensionful quantities are then measured in 
appropriate units of the lattice spacings. 

Note that Eq. [l]is a highly tuned action in the sense that two- and three-body interactions are absent. Radiative 
corrections giving rise to such operators may be avoided without fine-tuning of the lattice action through the use of 
a point-split four-body interaction, which disallows the formation of fermion loops on a fixed time slice. At infinite 
volume and for an attractive interaction, the binding energy E of the four-fermion system on the lattice may be 
related to the coupling by solutions to the integral equation: 



1 

2^g 



n 



27r 



Tier ^Pa 



(3) 



where ^p(m) = 1 + Ap/m, and Ap = 2sin^(p/2). This result is obtained by diagonalizing the four-fermion transfer 
matrix associated with Eq. [T] on the lattice at zero chemical potential, following the approach of . 

In the SU{A) symmetric limit (m^ = m for all a), the four-fermion coupling may also be related to the four-particle 
scattering length by evaluating the inverse scattering amplitude A'^^{p) for four-fermion scattering at vanishing 
external momentum p — \J mE, and equating the result to m/(47ra). This calculation may be performed on the lattice 



^ For a review of various other one-dimensional models for which the sign problem has been solved, see [22j and references within. 
^ Note that any choice of spatial boundary conditions are acceptable, however, the choice of open spatial boundary conditions simplifies 
the forthcoming analysis considerably. 
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following a prescription analogous to [2^, [25| , by summing the Feynman diagrams for four-particle scattering in a 
vacuum as a geometric series. Combining that result with Eq. |21 yields the relation: 

_ = 1 _ 1 . (4) 
Ana g 

The unitary limit (i.e., the limit a — oo) is obtained by tuning the coupling g to some critical value gc^ corresponding 
to a zero-energy four-body bound state (i.e., a solution to Eq. [3]with E = Q). At finite but asymptotically large 
positive scattering lengths, the binding energy of the system is given by i? ^ —\/{2ma?'). Note that negative values of 
the coupling are required when m/iA-Ko) exceeds l/(?c, and in conventional path-integral representations this usually 
implies a sign problem. As will be demonstrated in the proceeding sections, however, in the loop representation the 
sign problem is avoided for an expanded parameter regime governed by 5 > —1. 



III. LOOP REPRESENTATION 



The grand-canonical partition function for the lattice theory is given by 

Z{ii) = f [d^t][d^]e-5(M) ^ (5) 



where 

E 

neA 



(6) 



One may express Eq.[5]in terms of a hopping parameter expansion [26l |. by expanding e~^'^^^ in a power-series, followed 
by explicit integration over fermion fields. Note that at finite space-time volume, the expansion always truncates after 
a finite number of terms and is therefore well-defined, and furthermore converges (although the number of nonvanishing 
contributions will grow exponentially with volume). The contributions to the hopping parameter expansion from the 
various terms in the lattice action are given as follows: 



.ge^^'^ . (7) 



Upon integration over fermions, the surviving terms in the expansion are those in which all lattice sites n e A are 
occupied by a single power of ipt'^'u for each species a = {a, b, c, d}; multiple insertions of 4't'4'cr at the same site for 
a given species must vanish due to the ant i- commuting nature of Grassmann numbers. 

Nonvanishing configurations in the hopping parameter expansion may be specified diagrammatically in a unique 
way as a configuration of monomer and tetramer bonds associated with insertions (either zero or one) of either the 
nearest neighbor or four-fermion interaction at various locations of the lattice; for a given configuration, all remaining 
unoccupied lattice sites are saturated by insertions of the on-site interaction. Non-vanishing configurations under 
fermion integration satisfy the constraint that monomers form closed loops; for all such closed loop configurations, 
a tetramer may occupy time-like bonds in lieu of monomers, provided the bond was occupied by a monomer of 
every species. The contribution to the partition function from such configurations is simply the product of weights 
associated with each insertion of on-site, nearest neighbor and four-fermion operators, provided in Eq. [71 

Note that the hopping parameter expansion exhibits several convenient properties in one spatial dimension: 

• fermion loops formed by monomers of the same species are self-avoiding due to Pauli exclusion with one excep- 
tion: fermion loops may form on a fixed time-slice provided they only involve a single hop forward and back (or 
vice- versa); for every such loop, the configuration is weighted by a negative sign; 

• fermion loops may form by wrapping around the time direction and yield contributions to the partition function 
proportional to powers of e^^^; fermion loops may not wrap around the lattice in space, however, as a result of 
employing open spatial boundary conditions; 

• a fermion loop cannot wrap more than once around the lattice in time since closure of the loop would require 
the fermion path to cross itself. Multiple fermion loops are permitted provided each loop wraps around the time 
direction only once; the sign arising from every fermion loop of this kind cancels with the sign associated with 
anti-periodic temporal boundary conditions; 



• the number of sites that are saturated by insertions of the on-site interaction is even for every configuration 
provided Nt is even, therefore the sign associated with the weight for on-site interactions may be ignored. 



FIG. 1. Example of three class-equivalent closed-loop configurations belonging to the set C on an Nt x Ns = 6 x 6 lattice. 
The left-most configuration has <S = 8, = 16, Bt = 12 and T — A, whereas the center configuration has S = 10, Be = 14, 
Bt = 12 and = 3. The right-most configuration belongs to the subset C* of bubble-free configurations and has <S = 8, 
Bs = 8, Bt ~ 12 and J- = 0; this configuration has eleven domains d £ T>{c) (shaded) of which seven have = 1, three have 
C{d) = 2 and one has jC{d) = 3. 

With these considerations in mind, the partition function may be explicitly expressed as 

(l+gf^(^'^') , (8) 

where C is the set of all possible configurations involving self-avoiding, closed-loops which wrap around the time 
direction and space-like single-hop fermion bubbles. An example of two such configurations for a single species is 
shown in Fig. [T] (left, center). In Eq. [51 5(c) represents the total number of sites left unvisited by loops associated 
with c, Bs{c) and Br{c) are the total number of space-like and time-like bonds lying on the loops associated with c, 
and T{c) is the total number of space-like fermion bubble loops associated with c. Note that these quantities satisfy 
the relationship NrNg = Br + Bs + S; since Br (x Nr and Bs is always even, if Nr is even then so is S. Thus, I am 
justified in dropping the negative sign associated with the weight of the on-site interaction. 

Note that the quantity Br (c) /Nr is an integer, and may be regarded as a fermion charge associated with configuration 
c. Particularly, this quantity specifies the winding number of the configuration. Viewing a particular configuration 
c as a set of occupied bonds, Br{(^aCa) is defined as the number of time-like bonds in the intersection Ha-Ca- Such 
intersections may be associated with either four insertions of a monomer bond (one of each species) with a weight of 
unity or a single insertion of a tetramer bond with a weight of g, up to a common factor of e"^'''' which is accounted 
for separately. For every fixed set of configurations {c^}, I have summed over both possibilities, leading to an overall 
weight of 1 + 5 associated with the intersection on each time-like bond. 

In its present form, Eq. |8]has a sign problem associated with space-like vacuum bubbles on each time slice. Partic- 
ularly, the path-integral is weighted negatively whenever '^^J'ica) is odd. One may eliminate the sign problem by 
explicitly summing up all possible fermion bubble contributions. To see how this works, first note that the configu- 
rations c € C fall into classes identified by the path taken by closed-loops that wrap around the time direction; these 
classes are enumerated by the configurations c S C*, which involve no space-like fermion bubbles. An example of two 
configurations that belong to the same class are provided in Fig. [T] (left, center), along with their unique bubble- free 
class-equivalent configuration (right). 

Let me now carry out the sum over all configurations belonging to the same class. First, note that for every 
configuration c G C, we may identify a set of maximal one-dimensional space-like domains I?(c) associated with the 
sites left unvisited by closed loops, as demonstrated by the collection of shaded regions in Fig. [1] (right). Note that 
every configuration belonging to the same class as c shares a common set of domains V. Next, consider a single 
domain d G 2?(c) belonging to a fixed time slice and of linear spatial extent n = C{d), where < n < Ns- Since 
each domain is disjoint from all the others on every configuration belonging to the same class, we may assign an 
independent partition function Zn{m) to each of these domains. Specifically, I define Zn(jn) as a sum over the weights 
of all possible space-like fermion bubbles within a domain of length n. With such a definition, one may show that the 
total sum over all configurations belonging to the same class as c € C* will have a weight proportional to product of 
the partition functions Z£(c;)(™) associated with each of the domains d G T^^c). 

I evaluate Zn{m) explicitly for a given n by noting certain algebraic properties of the one-dimensional partition 
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FIG. 2. Diagrammatic representation for the recursion relation defined in Eq. [2]for 25(771). 



function. Particularly, the partition function satisfies the recursion relation: 



Zn{m) 



i(to) 1 + — 



1 



(2m) 



zo{m) 



1 



zi(to) 



1 

1 + — 

m 



(9) 



which may be represented diagrammatically as in Fig. [2] Note that the negative sign appearing in the second term 
of the recursion relation arises from the presence of a space-like fermion bubble loop. This recursion relation may be 
solved analytically and yields: 



z„(m) 



2«+Vl + 2A 




n+1 




n+1 



(10) 



one may easily confirm that this solution is positive for all to > and any n > 0. 0. 

Carrying out the summation over all configurations belonging to the same class as c € C*, I arrive at the final 
expression for the partition function: 



E 



n I n ^c{d){ma) 



2m„ 



(1 + 5) 



(11) 



where the sum over configurations is now constrained to the set of all possible self-avoiding fermion loops that wrap 
around the time direction unaccompanied by space-like fermion bubbles. This result is free of sign problems for any 
Too- > 0, jjLa and 5 > — 1, and is suitable for numerical simulation using a local loop- updating scheme or the "worm 
algorithm" (23], both at zero and finite temperatures. The canonical ensemble partition function associated with a 
fixed set of charges Q — {Qa, Qb, Qc, Qd} may be obtained via the Fourier transform of the grand-canonical partition 
function evaluated at imaginary chemical potential: 



Z{Q) 




(12) 



Under this projection, the configurations that contribute to the canonical partition function satisfy the fixed winding 
number constraint Qc — Br{ca)/Nr and belong to a subset of C*; an explicit expression for Z{Q) can be obtained 
from Eq. [TT] by simply replacing: 



,B^{c„)/N^ 



(13) 



IV. EXTERNAL POTENTIALS 

One may introduce a general space-time and species-dependent external potential v{x) in the continuum theory by 
replacing: /i — > /i -t- v{x) in Eq. [T] On the lattice , this corresponds to a modification of the discretization rules: 

(VV)n = /te6.(M-..)^^_^^ ^ (14) 

where = ^ (wn + v■a-eo)^ and n correspond to the half-lattice spacing locations n = n — ^eo for each n e A. As was 
the case with the chemical potential, an external potential is associated with the time- like bonds of the lattice and in a 
way may be viewed as a spacetime-dependent chemical potential. At first glance, this result may seem rather strange. 



^ This may also be seen by noting that Zn{m) = det [l — V'^/(2m)] , for a Laplacian defined on a spatial lattice of length n with open 
boundary conditions; the eigenvalues of the Laplacian are always real and positive. 
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since in the naive discretization, an external potential would be associated with sites, rather than time-like bonds. 
In fact, either choice is acceptable and should at most lead to different time discretization errors. Here, I choose the 
point-split discretization as it leads to a rather simple expression for the path-integral in the loop representation. 

In terms of a hopping parameter expansion, the introduction of an external potential suppresses fermion propagation 
in time at large potential energies. Hence, if and when fermion loops are present, they will most likely be localized 
about the minimum of the potential, to the extent that Pauli-exclusion will allow. In the loop representation, Eq. [11] 
is modified in an intuitive way, by replacing: 



where 



(15) 



(16) 



and the sum on n e c is over the half-lattice sites that lie on the time-like bonds of the configuration c. 

Provided that the external potential is confining, we may explicitly take the infinite volume limit without worrying 
about fermions wondering off to infinity. The reason we may do this is that the volume-dependence only appears in 
zc(d){m) for the domains d that touch the boundaries of the lattice. In any local updating scheme, only the ratios 



R^{m) = Zn±i{'rn) / Znim) matter; for the boundary domains, taking TV^ — ?> oo corresponds to taking the n 
of this ratio and yields: 




±1 



limit 



(17) 



where 



V. OBSERVABLES 

Short of explicitly introducing sources and sinks for the fields il}]^ and V'n, there are several observables which may 
easily be calculated using the above path-integral representation at finite temperature and density. Examples include 
derivatives of the partition function with respect to masses, potentials and the coupling, which yield expectation 
values for the kinetic energy operators (T^), potential energy operators (Kj), and interaction energy operator (/), 
respectively. More concretely, the expectation value of the kinetic energy operator is given by: 

iV.aiog(l/m.) - 91og(l/m.) + ''^^^"^ ] ^ 

aiog(l/m.) ^'^^ 

is a zero temperature and density contribution, which must be subtracted off. The expectation value of the potential 
and interaction operators are given by: 

\i ^ Y^aiogz(^ 1 

Nr ^ dlogWo-.fi 
n ' 

and 

1 9l0gZ(^) 1 _9.„ ^, 

Nr o log g Nr 1 + g 

respectively. In addition to these, charges and susceptibilities associated with each species are obtained from the 
expressions: 

Q.^^^logZ(,) = i-(^.(c.)), (22) 



and 



X.p = ^^a^f^log^^/^) = Wr [(^-(c-)^r(cp)) - {Br{c.)){Bricp))] ; (23) 
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these just measure the average winding numbers associated with the close loop configurations for each species, and 
their fluctuations. Combining these results, the total free-energy, F{fi) = — log Z (fi) / Nr , of the system in the grand- 
canonical ensemble is given by: 

F{p) = ^ (T, -t- K - fx,Q^) + I (24) 

(7 

at zero temperature. In the canonical ensemble, where the charges Qa are held fixed and the corresponding suscepti- 
bilities are zero (one should use an updating scheme designed to maintain such constraints), the expectation value of 
the total energy of the system is given by: 

£;(g) = ^(r, + K.) + / . (25) 

In the SU (4) symmetric limit, a quantity which is closely related to the interaction energy, known as the integrated 
contact density, is given by C = {mg)^I/g. This quantity appears in a set of universal relations, including the 
generalized Virial theorem: 

valid in the presence of a time- and species- independent harmonic potential = §(?^s — fij)^ centered about fis with 
spring constant k, and the adiabatic relation: 

' (27) 



d{-l/a) inm 

These relations may tested straight-forwardly in a numerical simulation using the loop representation for the partition 
function. 

Finally, density correlations may also be computed in the loop representation by considering functional derivatives 
of the partition function with respect to Vf^, the exponential of which couples to the density operator pf^ = -ip^e^ipn-eg- 
Generally, one may consider A^-point functions of the form: 

r^„...,,„(ni,...,nA,) = (^np,^,s,^ = l^n^^j log^(M) ; (28) 
if I define the local density field: 

0fi(c) = J2 ^^'i^ ' (29) 

then the number density of the system, n^- = rcr(n), is just given by the 1-point correlator: 

r,(n) = ((/)fi(c,)) . (30) 

Summing (n) over all space at an arbitrary time slice yields Eq. [22] for the expectation value of the charge of the 
system. Similarly the 2-point function (i.e., the density-density correlator) is given by: 

rai,<T2(ni,n2) = ((/'iii(c^i)(?!>fi2(CffJ) - (c^i )) (Ctra)) • (31) 

All other density correlators may be derived in terms of (c) from Eq. [511 in a straight- forward fashion. 



VI. MEAN-FIELD ANALYSIS 



The continuum limit for this theory corresponds to the simultaneous limits: a — > and — > 0, where a and rio- 
are measured in units of the lattice spacing, while keeping the dimensionless quantity rj = {ppa)~^ held fixed. Here, 
I investigate whether there exists such a regime in (g^ fj,„,m^) parameter space, where a continuum limit may be 
defined for any rj. I carry out the analysis within mean-field theory, deriving self-consistent equations for the number 
densities, and an upper bound /mf(m) on the free-energy density /(/i) = F{p.)/Ns, valid at finite temperature, for an 
arbitrary spin and mass imbalance, and for any coupling g > ^1. I numerically solve the self-consistent equations at 
zero temperature and study the phases of the model in the SU (4) limit of equal masses and chemical potentials for 
each species. 
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On general grounds, from the positivity of Eq. 1111 one may show that the free-energy density of the interacting 
four-component Fermi gas is bounded by j2&] : 



hiFif^) = /a (a*) + ^ {Sifi) - 5a (a^)), , (32) 



where f\{fJ-) is the free-energy density of some arbitrary reference system described by the action Sxil^-), the expectation 
value of the difference in actions is taken with respect to the reference system, and A is some variational parameter(s). 
The tightest bound on /(/i) is achieved by minimizing /a/fI/^) with respect to A. From Eq. 1111 the action associated 
with the partition function for the interacting Fermi gas in the loop representation is given by: 

5(m) = SpGit^a) - log(l + 9)Brin,c,) , (33) 

a 

where SpcifJ') is the action for a non-interacting single component Fermi gas at finite chemical potential. For the 
purpose of this calculation, I will choose as my reference action: 

S^ifi) = Spoif^-) - E ^Mca) ; (34) 

(7 a 

note from the form of Eq.[531 that S\{[i) = J2a SpGilJ-a + ^a), and therefore it follows that f\[y) — fpGii^a + ^a)- 
Evaluation of Eq. [32] with respect to the reference action yields 

fMF{y^) = E/FG(Mcr + Act) - log(l +5) ^'^^(Ma + A^) + ^ A^r JT-FG (A^ct + A^r) , (35) 
a a a 

where 



\ r r 1 / 

/fg(m) = !fg- - dp — log [l + e 



-A'^(log5j,(m)-Ai) 



(36) 



and 



fFG = -^l dp\og^p{m) ; (37) 



note that fpc is an irrelevant /^-independent shift in the free-energy density. To arrive at Eq. psp . I have used the 
relation: 

Br{n,C„)=J2l[^AiCa) : (38) 
n (J 

and the fact that an expectation value of the summand on the right-hand-side with respect to the reference system 
factorizes into a product of number densities associated with each species, each of which is independent of n at infinite 
volume. 

The mean-field free-energy density is minimized with respect to A when 

A, =log(l+5) [] 

a' ^(7 

where the free-gas number density is given by 

d 

npGip) = -^/j^g(m) ; (40) 
the mean-field density is in turn related to the free-gas number density by 

(T 

Specializing to the case of zero temperature^ the free-energy density for a single species becomes: 

fFGip) = /fg + SFGifJ') - A"t-fg(a') : (42) 



Here, I make use of the identity: lim^_>tx> log(l -I- x^)/ p = 9{x — 1) log(a;). 



10 
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FIG. 3. Left: number density phase diagram in the m — plane at ^ = 0, computed within mean-field theory. (I) and 
(IV) correspond to regions of vanishing number density, (III) corresponds to a region in which the lattice is fully saturated by 
fermions, and (II) is a region that is inaccessible by a mean-field treatment. A continuum limit can only exists in region (IV), 
where both the inverse number density (i.e., interparticle spacing) and scattering length can be taken large compared to the 
lattice spacing. Lines of separation between each of these regions are defined in the text. Right: number density as a function 
of chemical potential for a fixed fermion mass m — 5 and a — oo (i.e., starting from a point in region IV at ^ = 0). Dashed 
lines indicate extrema of the free-energy density, whereas the solid line indicates the solution to Eq. which minimized the 
free-energy density globally. The figure exhibits a first order phase transition from a region of small number density to a region 
in which the lattice is fully saturated by fermions. The continuum limit corresponds to the limit of zero density in lattice units, 
or /X ^ 0. 



where 

^FaifJ-) = - dp\ogip{m) , npcitJ-) = -ppifJ-) , (43) 

Jo 

are the corresponding energy and number densities defined in terms of the Fermi momentum: 



and ^0 = log (1 + 2/m) is the free-gas chemical potential at which the lattice becomes saturated by fermions. 

To simplify the analysis further, let me now consider the SU (4) limit in which the Fermi gas is unpolarized, and 
the mass of each species is equal. I begin by considering the number density phase diagram in the m — plane 
at /X = 0, represented by Fig. |3] (left). Within mean- field theory, this phase diagram is characterized by four regions, 
labeled (I), (II), (III) and (IV), separated by curves corresponding to the couplings: gi-ii — —1, gii-iii = 00, 
and gin -IV — e^'^^''^ — 1. Region (I) corresponds to repulsive interaction strengths ranging from gi-u < g < 0, 
and has vanishing number density within mean-field theory. Region (II) corresponds to a regime of strong repulsive 
interaction strengths ranging from —gii-m < g < gi-ii, and is inaccessible within a mean-field theory treatment 
due to the non-positivity of Eq. [TT] Region (III) corresponds to strong attractive interaction strengths ranging from 
giii-iv < g < gii-iii', in this regime, the number density is four, and corresponds to a lattice that is fully saturated 
by fermions. Finally region (IV) corresponds to attractive interaction strengths ranging from < g < gni-iv] in 
this regime the number density vanishes, allowing for the existence of continuum physics. 

Extending Fig. [3] (left) into a third dimension labeled by the chemical potential /i, the critical line defined by the 
coupling giii-iv becomes a critical surface separating a domain which is continuously connected to the zero density 
region (IV) from a domain which is continuously connected to the saturated region (III). In Fig. [3] (right), I plot the 
number density as a function of chemical potential for a fixed mass m — 5 and a — 00. Zero chemical potential on this 
plot corresponds to a point in region (IV) in Fig. [3] (left); as the chemical potential is increased, there exists a critical 
value at which the density changes discontinuously and one enters into the phase that is continuously connected to 
region (III) . Below this critical chemical potential, a region exists in which the density may be taken arbitrarily small 
and the scattering length arbitrarily large such that rj may take on a range of finite values, suggesting a continuum 
limit exists within the mean-field analysis. Note that the continuum limit places no restriction on the mass parameter; 
the mass m (or more generally, the reduced mass m = (X^o- l/'^o-) ^ in the mass-imbalanced case) may instead be 
regarded an anisotropy factor for the temporal and spatial lattice spacings. 
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VII. DISCUSSION 

I have demonstrated that the partition function for the four-component universal Fermi gas proposed in [TBj 
may be formulated on the lattice in a way that avoids the fermion sign problem at arbitrary chemical potentials, 
mass imbalances and sign of the four-body interaction strength. The formulation involved a path-integral over self- 
avoiding fermion loops, and is suitable for numerical simulation using standard local updating schemes. Furthermore, 
I have derived expressions for the expectation values of all relevant observables, including the kinetic, potential and 
interaction energies, moments of the charge distributions, as well as a general expression for N-point correlation 
functions involving the density operator. These observables i nay be directly related to interesting nonperturbative 
quantities in the BEC and cross-over regimes, as discussed in [15l|; full-scale numerical simulations to compute these 
quantities will be the subject of a future study. 

Positivity of the path-integral in the loop representation permits a mean-field analysis of the lattice theory. In the 
SU{A) limit, I have demonstrated within this analysis that there exists a phase in which both 1/a and n may be 
taken to zero (in lattice units) while achieving a range of finite values for 77, thus allowing one to define a continuum 
limit for the lattice theory. It would be interesting to extend this mean-field analysis to the population- and mass- 
imbalanced interacting Fermi gas, for which the path-integral representation retains its positivity property; ultimately, 
a nonperturbative study of the {g, m^) phase diagram is desirable as well. 

The formalism presented here may be generalized to include two- and three-body contact interactions in a straight- 
forward manner. Extending the approach to higher dimensions, while straight-forward, yields an action with a sign 
problem, however. The reason is that in higher dimensions, a single fermion loop may wrap more than once around 
the lattice in the time direction without crossing itself. Generally a negative sign accompanies every fermion loop, 
and an additional negative sign accompanies every instance that the fermion wraps around the lattice in time, due to 
the choice of anti-periodic boundary conditions. A relative negative sign in the weight will therefore always persist 
whenever the winding number minus the total number of fermion loops wrapping in time is odd. 

Finally, it was observed that in one spatial dimension, one may regard a nonrelativistic Fermi gas as a Bose gas with 
an impenetrable hard core [2^ . In light of this, it may seem unsurprising that the Fermi system may be formulated 
in a way that avoids the sign problem if the corresponding bosonic system is also free of sign problems. It is known, 
however, that the conventional path-integral formulation for bosons does in fact have a sign problem at finite density. 
Interestingly, the sign problem for bosons at finite density may be avoided by resorting to techniques similar to those 
discussed here [SQ]. 



ACKNOWLEDGMENTS 



I would like to thank J.-W. Chen, D. B. Kaplan, Y. Nishida, D. T. Son, and H. Suzuki for interesting and helpful 
discussions. I am grateful to the national Institute for Nuclear Theory, where portions of this work took place. 



S. Giorgini, L. P. Pitaevskii, and S. Stringari, [Rev. Mod. Phys., 80, 1215 (2008) 

G. Baker, in Recent progress in many-body theories: the proceedings of the 10th international conference, Seattle, USA, 
September 10-15, 1999, Vol. 3 (World Scient ific Pub Co Inc. 2000) p. 15. 
S. Gandolfi, K. E. Schmidt, and J. Carlson, 
M. M. Forbes, S. Gandolfi, and A. Gezerlis 



, _ _ , , , hys. Rev. Lett.. 106, 235303 (2011) 

J. Carlson, S. Gandolfi, K. E. Schmidt and S. Zhang, |Phys. Rev. A, 84, 061602 (2011) 



hys. Rev. A, 83. 041601 (2011' 



Phys. Rev. Lett.. 106, 205302 (2011) 
100, 150403 (2008)|. 



J. E. Drut, T. A. Lahde, and T. Ten, 

J. Carlson and S. Reddy, Phys. Rev. Lett. 

K. B. Gubbels and H. T. C. Stoof, ArXiv e-prints (2012), arXTv: 12 05.0568 [cond-mat.quan t-gasJ| 
G.-B. Jo. Y.-R. Lee. J.-H . Choi, C. A. Christensen, T. H. KiniT^J. H. Thywissen, D. ErPritchard, 



Science, 325, 1521 (2009) 



and W. Ketterle, 



http://www.sciencemag.or g/content/325/5947/1521.fulI.pdf 
Partridge. Y. A. Liao. and R. G. H ulet. International Journal of Modern Physics B, 23, 3195 (2009) 
Nuclear Physics A. 210. 157 (19731, 



W. Li, G. B_ 
V. Efimov 

E. Braaten and H.-W. Hammer 
I. Bloch, J. Dalibard 
Y. Nishida and S. Tan 



ISSN 0375-9474. 



Physics Reports. 428. 259 (2006\ ISSN 0370-1573 



Nishida and D 
Tan, 
Tan, 
Tan. 



and W. Zwerger. [Rev. Mod. Phys.. 80. 885 (2008) 
Phys. Rev. Lett.. 101. 170401 (2008)[. 



Annals of Physics. 323. 2952 (2008) 



T. Son, phys. Rev.. A82 043606 (2^151 7 arXiv: 0908 . 2 159 [cond-mat .quant-gas]] 



Annals of Physics. 323. 2971 (2008^ 



Annals of Physics, 323, 2987 (2008) 



cond-mat/0505200 

cond-mat/0508320J 

arXiv:arXiv:0803.0841 [cond-mat. stat-mech] [ 



12 



E. 

E. 

Y. 

H. 

J.- 

M 

D. 

C. 

N. 

J. 

M 

M 



Rev. Lett.. 100, 205301 (2008^ 



Braaten, D. Kang, and L. Platter 

Braaten and L. Platter, |Ph^^^^^_^_^^_^^^^^^^^ 

Nishida an d D. T. Son. [Phvs. Rev. D. 76. 086 004 (2007)^ 

G. Evertz, [Advances in Physics, 52, 1 (2003? " 

W. Chen and D. B. Kaplan 

, G. Endres, D. B. Kaplan, J.-W. Lee 

B. Kaplan, M. J. Savage, and M. B. Wise, [Nuclear Physics B 



Phvs. Rev. A. 78. 053606 (2008)[ , ' arXiv:0806.2277 [cond-mat .other]! 



Phys.Rev.Lett., 92, 257002 (2004) 



http://www.tand fonline.com /doi/pdf/ 10 . 1 080 /GOO 1873021000049195 
arXiv:hep-lat/0308016 [hep-lat] 
T2012), arXiv: 1203 .3169 [hep-lat 



and A. N. Nicholson, 



478, 629 (1996)1 arXiv:nucl-th/96050C^ 



Itzykson and J. M. Drouffe, Statistical Field Theory (Cambridge Univers ity Press, Cambridge, 1989) 



Prokof ev and B. Svistunov, [Physical Revi ew Letters, 87, 160601 (2001)[ arXiv:con d-mat/0103146, 
Jensen, Acta Mathomatica, 30, 175 (1906), ISSN 0001-5962 10.1007/BF02418571. 



Girardeau, Journal of Mathematical Physics, 1, 516 (1960) 
G. Endres, |phys. Rev. D, 75, 065012 (2007)[ . 



